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Abstract 

We propose explicit recipes to construct the euclidean Green functions of gauge- 
invariant charged, monopole and dyon fields in four-dimensional gauge theories 
whose phase diagram contains phases with deconfined electric and/or magnetic 
charges. In theories with only either abelian electric or magnetic charges, our 
construction is an euclidean version of Dirac's original proposal, the magnetic 
dual of his proposal, respectively. Rigorous mathematical control is achieved for 
a class of abelian lattice theories. In theories where electric and magnetic charges 
coexist, our construction of Green functions of electrically or magnetically charged 
fields involves taking an average over Mandelstam strings or the dual magnetic 
flux tubes, in accordance with Dirac's flux quantization condition. We apply our 
construction to 't Hooft-Polyakov monopoles and Julia-Zee dyons. Connections 
between our construction and the semiclassical approach are discussed. 
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1. Introduction 



In this paper, we study a variety of gauge theories in four dimensions with the 
following features: their phase diagrams contain phases with deconfined electric or 
magnetic charges; their particle spectra thus contain electrically or magnetically 
charged particles or dyons. They may exhibit phase transitions characterized by 
condensation of charged particles or magnetic monopoles. Some of them exhibit 
duality symmetries. 

One is interested in studying various aspects of the phase diagram and of the dy- 
namics in particular phases of such theories. Obviously, one would like to study 
these theories analytically. However, explicit analytical results are often only avail- 
able for physically rather unrealistic theories with super symmetries. In order to 
study more realistic theories, one therefore often resorts to numerical investigations 
of lattice approximations of such theories. Usually, such investigations are based 
on a euclidean (imaginary-time) formulation of quantum field theory obtained 
from a real-time formulation by Wick rotation. In the lattice approximation, one 
replaces Euclidean space-time by a (finite, but arbitrarily large) lattice. 
In order to study the mass spectrum of a quantum field theory, one considers 
euclidean Green functions of gauge-invariant (physical) fields of the theory that 
couple a vacuum (ground) state to some one-particle state. Masses can be calcu- 
lated from exponential decay rates of certain two-point euclidean Green functions. 
Signals for a phase transition, e.g., one between a Coulomb- and a confining phase, 
can be detected in an analysis of asymptotic behaviour of suitable euclidean Green 
functions. For example, the transition from a Coulomb- to a confining phase in 
an abelian gauge theory is reflected in the appearance of long-range order in the 
two-point euclidean Green functions of magnetic monopoles. 
On a more foundational level, we are longing for a complete description of quan- 
tum field theories in terms of the (euclidean) Green functions of gauge-invariant 
interpolating fields. For example, to describe the deconfined phase of an abelian 
gauge theory, we would like to construct Green functions of charged fields. 
For QED-like theories without dynamical magnetic monopoles, a proposal for a 
gauge-invariant charged field has been made, many years ago, by Dirac [1]. For 
lattice theories of this type, Dirac's proposal has been studied carefully e.g. in 
[2,3]. For the convenience of the reader, the main results of this analysis are 
summarized in sect. 3. 

There are, however, plenty of physically interesting theories with deconfined , elec- 
trically charged particles and dynamical magnetic monopoles. The problem of con- 
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structing gauge-invariant electrically or magnetically charged interpolating fields 
for such theories has not been solved in adequate generality and is quite non- 
trivial. The difficulties encountered in studying this problem are a consequence of 
the Dirac quantization condition. We report partial results towards a solution of 
this problem in sect. 4. 

The main purpose of this paper is to describe constructions of gauge-invariant elec- 
trically or magnetically charged interpolating fields and of dyon fields in (lattice) 
gauge theories with dynamical electric charges and dynamical magnetic monopoles 
and to exhibit duality transformations converting electrically into magnetically 
charged fields (and vice versa). Our work is primarily kinematical : we propose 
fairly explicit recipes for how to construct the euclidean Green functions of such 
fields. But we do not engage in any mathematically careful analysis of the prop- 
erties of these Green functions. Instead, we gather evidence supporting various 
conjectures on their behaviour. Part of this evidence comes from previous, math- 
ematically precise work on lattice gauge theory, another part is based on more 
heuristic arguments and conventional wisdom. 

All in all, we believe we arrive at a fairly consistent picture of how electrically or 
magnetically charged and dyon fields should be constructed. 

We also review duality properties of some gauge theories and analyze how duality 
transformations act on charged and dyon fields. This part of our analysis makes 
contact with issues that have been quite topical, during the past few years; (see 
e.g. [4]). 

Finally, we outline a formal construction of Green functions of monopole- and 
dyon fields in the continuum SU(2) Georgi-Glashow model and indicate how our 
construction is related to the (semi-) classical analyses of 't Hooft and Polyakov 
[5] and of Julia and Zee [6]. 

Next, we present brief summaries of the different sections of this paper. 
In section 2, we introduce the gauge theories studied in this paper. We consider 
three classes of models (A, B, and C). The models in class A are non-compact 
abelian lattice gauge theories with electrically charged matter fields, but with- 
out dynamical magnetic monopoles, i.e., models of lattice QED. The models in 
class B are compact abelian lattice gauge theories with electrically charged matter 
fields. As originally pointed out by Polyakov [7] , they describe dynamical magnetic 
monopoles coexisting with electrically charged particles. The models in class C 
are related to (lattice approximations or formal continuum limits of) the Georgi- 
Glashow model. We define our models in terms of euclidean action functionals. 
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We introduce some key notation and review some facts on the phase diagrams of 
these models. 

In section 3, we construct gauge-invariant charged fields for models of class A, in 
accordance with Dirac's proposal and with the result of [2,3]. We also construct 
monopole fields in compact, pure abelian lattice gauge theories. These theories are 
dual, in the sense of Kramers- Wannier duality [8], to some models of class A, and 
this suggests a dual variant of Dirac's proposal as the right definition of monopole 
fields. 

We start section 3 with a short recapitolation of Osterwalder-Schrader reconstruc- 
tion and of an analysis of superselection sectors based on euclidean Green func- 
tions of charged fields and monopole fields. We then recall (an euclidean version 
of) Dirac's proposal [1] for gauge-invariant charged fields, 

E) = $ x e tiA ' E{x)) (1.1) 

where <& x is a charged matter field, and A is a (non-compact, i.e., real-valued) 
abelian gauge potential; furthermore, E(x) is a c-number one-form related to 
the electrostatic Coulomb field of a point charge. (Thus, E(x) has a source at the 
point x whose charge is the same as the charge of $). We study the "infra-particle 
nature" of charged particles in such theories. 

Our construction of monopole fields and our analysys of their Green functions in 
compact abelian gauge theories without matter fields follows from the results on 
electrically charged field by using Kramers- Wannier duality; (this is made explicit 
in sect. 3.3). 

In sect. 4, we study electrically charged particles and dynamical magnetic monopoles 
in compact abelian lattice gauge theories with matter fields. We start by explain- 
ing the origin of Dirac's quantization condition, 

q e -q m = 2ym (1.2) 

n = 0, ±1, ±2, .... We then construct gauge-invariant electrically charged fields as 
averages of charged fields multiplied by Mandelstam string operators (exponential 
of the gauge field integrated along a path, called Mandelstam string [9]), the av- 
erage being taken over a suitable space of Mandelstam strings. Among the more 
subtle points appearing in this paper is the one to come up with a good defini- 
tion of "averages over Mandelstam strings". In an analysis of Green functions 
of these fields based on successively integrating out the large-frequency (short- 
distance) modes of the fields, one observes that the large-distance effective theory 
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becomes increasingly similar to the one describing a model in class A. In partic- 
ular, the functional integral expressions for these Green functions resemble the 
ones for the Green functions of charged fields constructed according to Dirac's 
proposal, in models of class A. Suitable averaging over Mandelstam strings yields 
operators which, under renormalization, approach ones related to (1.1). The mag- 
netic monopoles get suppressed, and Dirac's flux quantization condition becomes 
irrelevant. 

We then proceed to define monopole fields, or, more precisely, euclidean Green 
functions of such fields, in terms of certain averages over 't Hooft disorder operators 
[10]. We also present a heuristic analysis of properties of euclidean Green functions 
of electrically charged- and monopole fields and compare our results with these in 
section 3. 

We conclude section 4 with an analysis of dyons and of an SL{2, Z) duality group 
in a compact abelian lattice gauge theory with a topological term, related to the 
instanton number, in the action. Our analysis is based on previous work of Cardy 
and Rabinovici [11]. 

In section 5, we outline a construction of physical interpolating fields for 't Hooft- 
Polyakov monopoles and Julia-Zee dyons in the Georgi-Glashow model on the 
lattice and in the formal continuum limit of this model. The role of 't Hooft's 
7*2 monopoles and the corresponding disorder operators in our construction in 
explained in detail. 

Some technical points are relegated to two short appendices. 
2. The models 

The models we consider in this paper are lattice gauge theories. ( Some comments 
on a continuum gauge theory are added in the last section, for the Georgi-Glashow 
model.) 

Our (euclidean space-time) lattice is Z| , where the subscript | indicates that the 
coordinates of the sites are half-integers. We will also have to consider sublattices 
(more precisely cell subcomplexes) of Z\ . 

2 

Lattice fields can be defined as follows: 

A scalar field $ is a map from the sites, i, of the lattice to a normed vector space 
V H . 

A fermion field \1/ is an anticommuting map from the sites to the orthonormal 
frames of a vector space V s <E> Vf, the fermion space, where V s , the spin space, 
carries a representation of the Dirac-Clifford algebra. 
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A gauge field g is a map from the links < ij > of the lattice to a group G, the 
gauge group. 

If W is an additive abelian group and k is a positive integer, one can define k— 
forms with values in W as maps, F k , from oriented fc-dimensional cells, c k of the 
lattice to W satisfying F(— c k ) = — F(cfc), where denotes the cell obtained from 
Cfc by reversing the orientation. 
We denote by d the lattice exterior differential 

dF(c k+1 ) = F (°^ 

c k edc k+1 

and by * the Hodge-star: 

Let c\_ k denote the cell in the dual lattice, Z 4 , dual to c k . Then 

*F(ct_ k ) = F(c k ) 

We also introduce the co differential 5 = (— ) d ( fc + 1 )*d* and the Laplacian A = 
d8 + 8d. If A is a sublattice we denote the restriction of lattice operators to forms 
over A by a subindex A, e.g. d\ instead of d. 

If W is a Hilbert space with scalar product (•, •) one can define a scalar product 
among k— forms F and F' by 

(F,F') = J2(F(ck),F'(c k )) (2.1) 

Cfc 

The restriction of the scalar product (2.1) to forms defined on a sublattice A is 
denoted by (■, -)a- We define the £ 2 -norm of F by 

\\F\\ = yf(F,F) 

The vacuum functional of a gauge theory with gauge field g and with matter fields 
$ is given in terms of formal integration "measures" 

dfji(g, = i J] ^ <i7> n rf$ l ^,rfvl/ l e- s ^'*' ,|r ' >1 ') (2.2) 

<ij> i 

where dg < ij > is the Haar measure on G, <i$^ is the Lebesgue measure on Vh, d^id^i 
denotes Berezin integration, S is the total action, and Z is a normalization factor, 
the partition function. [Mathematically, the measure (2.2) is first defined for a 
finite lattice A C Z| , with suitable boundary conditions at <9A; subsequently, one 
takes the limit A / Z\\. 

2 
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Expectation values w.r.t. the measure (2.2) are denoted by (•). 
2.1 Actions 

We consider the following classes of models. 

Class A: Non-compact abelian gauge theories with matter fields . 

As an example we discuss the abelian Higgs model . 

The fields of this model are a real gauge field, A, and a complex scalar field, $. 
The action is given by 

S(A, d>) = S (A) + S 1 (A, d>) + S gf (A) 
S (A) = ^\\dA\\ 2 , S 1 (A, *)=£2><- e M <->^| 2 + A^(|$,| 2 - l) 2 

<ij> i 

(2.3) 

Sgf(A) is a gauge fixing term, e.g. 



S 9f (A) = 




From now on, this gauge fixing term will usually not be written explicitly in our 
formulas, anymore. In the limit A /* oo this model reduces to the Stiickelberg 
model. For simplicity, we only discuss the model at A = oo, in this paper. It will 
be referred to as model A. 

Remark 2.1 A second interesting model in class A is spinor QED on the lattice; 
(see [2] for a rather detailed discussion). The fields of QED are the real gauge field 
A and a four-component fermion field ^ = {^ a , a = 1, ...4}. The action is given 
by 

S(A, V, *) = S (A) + S 1 (A, V, *) + S gf (A) 

<ij> i 

where r < jj > are matrices on Vp given by 

r<v> = i ± 7/i 

if < ij > is directed in the ±fi direction, fj, = 1, d, where 7 M are Euclidean Dirac 
matrices. 

Remark 2.2 In these models, we can omit gauge fixing by working directly with a 
measure defined on gauge equivalence classes 
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[A] = {A':A'-A = d£}, 



(2.5) 



where £ is a real scalar field [12]. 

Class B: Compact abelian gauge theories with matter fields . 

As an example we analyze the compact abelian Higgs model (model B). 

In the "Villain formulation", the basic fields are a i7(l)-valued 0-form (p, a £7(1)- 

valued 1-form 6>, a 27rZ-valued 1-form £ and a 2-7rZ-valued 2-form n. The action 

is given by 

S(6,n,<p,e) = S (6,n) + S^O^J) 

s (e,n) = ^\\de + n\\ 2 , Sl (e,<fi,e) = ^\\d<f + qd + e\\ 2 , (2.6) 

where q is the charge of the matter field; it is an integer. 
Alternatively, we can use the "Wilson formulation" : 

$>(*) = § E^" 003 ' S 1 (9,<p) = % J^(l-cos(d<p + q e) <tJ> ). 

P <ij> 

(2.7) 

Class C: Non abelian gauge theories coupled to matter fields breaking the gauge 

group to a Cartan subgroup containing U(l) . 

As an example we analyze the Georgi-Glashow model (model C). 

Its basic fields are an SU(2)-valued gauge field g and a Higgs field <E>, of unit 

length, in the adjoint representation of SU(2). The action is given by 

S(g,Q) = S (g) + S 1 (g,$) 
S (g) = P^i 1 ~ Xi9dp)) , S 1 ( S ,*) = 0h 1 Z (^U H ( g<ij> )<S>) (2.8) 

P <ij> 

where x denotes the character of the fundamental representation, and Uh(-) de- 
notes the adjoint representation. 

2.2 Phase diagrams 

For small k and model A has a Coulomb phase with a massless photon 

[13,12,14]. Furthermore, it is known that, for small (3 and suficiently large k, it 
has a superconducting Higgs phase [15]. 
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Model B, with q = 1, has a confining / Higgs phase, for (3 small or k large . 
Furthermore one expects a massless phase for f3 and k small [16]. 
For q 7^ 1, a phase transition line is expected to separate the confining and the 
Higgs phases, furthermore, for q sufficiently large, a Coulomb phase is known to 
appear [13,12] for intermediate values of for arbirary large k. If we add a 
topological 6-term, 



where the wedge product on the lattice can be defined as in [17], to the action of 
model B, in the limit k /* oo, we recover the Z q model discussed by Cardy and 
Rabinovici, which is expected [11] to exhibit Higgs, Coulomb, confinement and 
"oblique confinement" phases , related to an (approximate) symmetry under an 
SL{2, Z) duality group. 

Model C has been rigorously shown to have a Coulomb phase for f3 h = oo and /? 
large. The Coulomb phase is expected [18] to extend into a region of large values 
of /3h and f3. A confining phase is also known to exist [15], for small values of (3. 
Charged particles and monopoles are expected to exist in the Coulomb phase of 
these models as (infra-) particle excitations. This is the main issue discussed in 
this paper. 

For models in class C, magnetic monopoles are expected to exist as massive par- 
ticles also in the continuum limit (provided it exists). These are the celebrated 
at't Hooft-Polyakov monopoles. In the models of class B, monopoles are a lattice 
version of the (singular) Dirac monopoles. 

3. Charges or monopoles 

To motivate our constructions of charged and monopole fields for the models 
introduced in the previous section, we outline a general strategy exploited in 
[19,20,21,2] (see also [22]) to construct charged (and soliton) fields and superse- 
lection sectors directly from correlation functions defined through euclidean func- 
tional integrals. We discuss this construction on the lattice; but similar ideas have 
been applied to continuum models, although only formally. 

3.1 Reconstruction of charged and soliton quantum fields from euclidean Green 
functions 

Let O(C) denote a euclidean "observable", i.e., a neutral, gauge-invariant, local 
function of the basic fields with support on a compact connected set of cells C, 




(2.9) 
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such as a Wilson loop, a " string field" , etc. 

We consider expectation values of such euclidean observables 

G m {C U ...,C m ) = (0{C 1 )...0{C m )). 

Let r denote reflection in the time-zero plane, / a complex-valued function and / 
its complex conjugate. We define the Osterwalder-Schrader (O.S.) involution ©os 
as an antilinear map satisfying 

©Os/($i) = f($ri),®Osf(g<ij>) = f(g<rtrj>), 

©os** = *ri7o, ©os 1 ^ = 7o*™ (3.1) 

Let F + denote the set of linear combinations of euclidean observables 0(C) with 
C supported in the positive-time lattice. If, for every F G F+, 

((QosF)F) > 

then the correlation functions {G m } are said to be O.S. positive. 

Theorem 3.1 (O.S. reconstruction). If the correlation functions {G m (Ci, C m )} 

satisfy 

i) lattice-translation invariance, and 

ii) O.S. positivity 

then one can reconstruct from {Gm} 

a) a separable Hilbert space, Ho of physical states, 

b) a vector of unit norm, O G Ho, the vacuum, 

c) a self-adjoint transfer matrix T, with unit norm, and unitary space translation 
operators U^, n = 1, ...d — 1, acting on Ho and leaving Q invariant. We define 
T(t) = T*, t G Z + , U(a) = Ut=i U^,S= (a l7 a d _0 G Z^" 1 . 

If, moreover, the correlation functions satisfy 

iii) cluster properties 
then 

d) Q is the unique vector in Ho invariant under T and U^. 

From the explicit proof of the theorem it follows that there is a set of vec- 
tors {\Ci, ...,C m >} in Ho, with Cj contained in the positive time lattice, j = 
1, m,such that the set of linear combinations, denoted F+, is dense in Ho- On 
these vectors, the scalar product is defined by 
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(C 1 ,...,C m \C[,...,C' n } = (Oos[0(C 1 )...0(C m )]0(C[)...0(C' n )} (3.2) 

Field operators 0(C\, ...,Cg), with Cj contained in the strip {x : x° E (0,t),t G 
Z_|_}, can be defined on T{i)F+ by setting 

6(C U C e )T(t)\C[, C' m >= |d, C £ , (C[) t , (C'J t > (3.3) 

where (-) t denotes translation by t in the time direction. The operators 0(C) are 
the quantum mechanical operators reconstructed from the euclidean observables 
0(C). The algebra, A, generated by the operators 0{C\, Cg)T(t) defined above 
is called " lattice observable algebra" . 

If the space of the physical states of the model contains charged or soliton states, 
such states are not contained in the Hilbert space constructed above. 
The construction underlying Theorem 3.1 reproduces only the vacuum sector of 
such models. Suppose, however, that the Hilbert space, TL, of physical states of the 
model can be decomposed into orthogonal sectors Ti q invariant under T, U ^ and 
the lattice observable algebra A, i.e., 

T~C = ®ql~tq- 

Here Hq is defined to be the subspace given by Afl. It is called vacuum sector . 
A sector H q _L Ho is called a "charged sector" if there are no lattice translation- 
invariant vectors in TC q . 

To construct the charged sectors, suppose that we can find an enlarged set of cor- 
relation functions {G n>m } obtained by taking expectation values involving, besides 
euclidean observables, "charged" order- and disorder-fields, denoted by C(T) and 
by T>(T), respectively, with support in a connected (but, in general, non-compact) 
set of cells V. Correlation functions with non vanishing total "charge" q are ob- 
tained from those of vanishing total charge, by removing the support of the charge 
of a "compensating" field of charge — q to infinity. 

Assume that the correlation functions {G n>rn (Ci, C n , Ti, T m )} still satisfy 
the hypotheses of the O.S. reconstruction theorem. 

Denote by 7Y,T, Q,,0(C),C(T),T>(T) the Hilbert space, the transfer matrix, 
the translation operators, the vacuum and the quantum fields obtained via O.S. 
reconstruction from the correlation functions {G n>rn \. If all the correlation func- 
tions of non- vanishing total "charge" vanish, and clustering holds, then the Hilbert 
space H contains charged sectors, because the scalar product between two charged 
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states of unequal charge, defined in analogy with (3.2), vanishes; (for more precise 
definitions and proofs see [ ]). 

The quantum fields C(T), V(T) map the vacuum sector to a charged sector. 

Remark 3.1 Some of the quantum fields might couple the vacuum to one-particle 
states, or infra-particle states. This can be seen as follows: Provided that T > 0, 
the mass operator can be defined by 

M = -ln(T\{H i0) eCfi}) (3.4) 

where Ti.^ denotes the fibre of H. of zero total momentum, (i.e. a generalized 
vector |$ > belongs to iff >= |^ >). Then we have the following 

result. 

Theorem 3.2 A field operator A acting on Ti couples the vacuum O to a stable 
massive one-particle state iff 



e -m(A)t 

< AQ,T(t)U(a)An > - < Afl,Q >< Q,AQ > as t / oo (3.5) 

t 2 

with m(A) > 0, for any a G Z d_1 . This decay law is called Ornstein-Zernike 
decay. (For infra-particles the exponent in the denominator on the r.h.s. of (3.5) 
has a positive small correction.) 

For many lattice models, the large t behaviour involved in (3.5) can be analysed 
in terms of expansions methods [23,19,21]. 

3.2 Charged Gelds in non-compact abelian models 

In this section we recapitulate the basic steps of the construction of charged fields 
in the models of class A [2,3], following the scheme outlined above. We do this 
for the sake of completeness and in order to elucidate the difficulties arising in 
attempts to extend our construction to models in class B. 

Let E be a real-valued 1-form with support on an infinite, connected sublattice 
of the time-zero hyperplane,Ao, of the lattice Z 4 . 
Furthermore we assume that 



5E = d (3.6) 

where 



,g \ (1 if « = 

1 otherwise 
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and 

E<ij> ~ d(< ij >, 0)~ 2 , as d(< ij >, 0) / oo (3.7) 

Here is the origin in Z 4 and d(< ij >, 0) denotes the euclidean distance between 
< ij > and 0. As an example, one may consider 

E = d Ao AX> (3.8) 

We denote by E(x) the 1-form E translated by x; E(x) describes the electrostatic 
Coulomb field surrounding a source of charge 1 located at x. We define the charged 
fields of the Higgs model by 

ft E) = ($ x y e iq ( A ' E(x) \ ${x, —ft E) = ($ x ) q ' e -^,E(x)) ^ q £ N 
(See fig.l) Typical observables, O, are Wilson loops 

Q{aC) = \\ e iA<i i> a (3.10) 
<ij>ec 

where C is a loop, a G R/{0}, and "string fields" 

0{C xy ) = * x H e iA <^>% (3.11) 
<ij>ec xy 

where C xy is a path form x to y. 
Correlation functions G^ m are defined by 

n m 

Gl m (x iqu x n , q n , C u ...C m ) = (J] ft, E) J] 0(C,)) (3.12) 

i=i 3=1 

For later purposes it is convenient to define a generalization of (3.12): For two 
different 1-form E, E' satisying (3.6), we set 

G n+ m, r (xi, ft, E, x n , q n E, x[, q[, E', x' m , q' m , E' , C\, ...C r ) = 

n m r 

= qi s) n ^ ft', so n 0(0) (3.13) 

i=l j = l ^=1 

The correlation functions G n ^ m can be expressed as sums over configurations of 
"electric currents" , by first integrating over the matter fields and then integrating 
over A. For example 
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V , z(o)e~% (p+ e > a ~ (p+ e )) 

where p is an integer valued 1-current, z(p) is a certain statistical weight deter- 
mined by the action of the model, and 

n 

E = Y J <liE{xi) 

i=i 

Since, in the denominator of (3.18), dp — 0, the support of the currents p is given by 
a set of loops. Hence the denominator can be interpreted as the partition function 
of a gas of current loops interacting via the four-dimensional (lattice) Coulomb 
potential A -1 . In the numerator of (3.18), open currents appear, with sources 
at {xi}, besides current loops. At the sources, the electric currents spread out in 
fixed time planes, as described by E{xi). The currents p can be interpreted as the 
Euclidean worldlines of charged particles; currents supported on loops correspond 
to worldlines of virtual particle-antiparticle pairs, open currents with connected 
support correspond to the worldlines of particles created at one end of the line 
and annihilated at the other one. 

Theorem 3.3 The correlation functions G E m defined in (3.12) are lattice transla- 
tion invariant and O.S. positive. Furthermore, for (5 large enough or for strictly 
positive j3 and k small enough, clustering holds and all correlation functions with 
non-zero total charge vanish. 

Idea of proof Invariance under lattice translations follows from the existence of the 
thermodynamic limit of the measure corresponding to (2.3) derived by correlation 
inequalities [15,16]. O.S. positivity of those measures can be proved as in [15], and 
O.S. positivity of charged correlations follows from the fact that E(x) is localized 
in a fixed-time plane, so that e.g. 

G E (rx,-q,y,q) = (G s ((<&.) V^'^) (* y )*eW A > E <»») 

for x°, y° > 0. (There is a slight subtlety cocerning the choice of boundary condi- 
tions for correlation functions of charged fields in a bounded space-time volume. It 
can be dealt with in a way similar to that explained in [21]; see also section 4.2) 
Cluster properties are a consequence of (generalizations of) the bounds stated in 
the next theorem, for details see [2,3]. 

Theorem 3.4 In model A, for f3 large enough or k sufficiently small, and \x — y\ 
large enough, 
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G 2)0 (x, -1, E, y, 1, £■') < exp{- Cl (/3, £)} 

< exp {- (| - c 2 (/3, «)) ((£ - E dip )^-\E - E dip j) }c 3 (p, k) ^ (3.15) 
and for x = ry 

exp + pmin^A-^E + /w))}c(K) |x - y| < G 2i0 (x, -1, £, y, 1, £') (3.16) 

where Cj(/3, k) tends to exponentially, as /3 /* oo, and, for fixed /?, Ci((3,K) ~ 
O(k), as « | 0, z = 1, 2, 3, and c(k) < 1, c(«) ~ 0(k), as k | ; Pmin is a current 
of minimal length and flux 1 connecting x to y, E = E(x) — E'(y) and 

E dip = dA- 1 (6 x -6 v ) (3.17) 

The upper bound in Theorem 3.2 shows that charged (infra-)particles in the Higgs 
model have strictly positive mass in the range of coupling constants indicated in the 
theorem. The lower bound proves that, in the same range of coupling constants, 
their mass is finite. 

The method of proof is based on a combination of a Peierls- and renormalization 
group argument, following [12,24]. For k = oo the lower bound in (3.16) can be 
obtained more easily from Jensen's inequality, using representation (3.14). 

Remark 3.2 It can be shown, following [25], that, for f3 sufficiently small and k 
sufficiently large, clustering holds, and 

G 2 , (x, 1, E, y, -1, E') > e -omE,(A+o m -i)E) > CQnst (3 Jg) 

uniformly in \x — y\. Hence correlation functions of non zero total charge do not 
vanish. This is a manifestation of charged-particle condensation typical for the 
Higgs phase. 

From the O.S. reconstruction theorem we obtain a Hilbert space, denoted by H(E), 
a transfer matrix T E and a dense set of vectors: 

\x 1 ,q 1 ,...,x n ,q n ,C 1 ,...,C m > E , q t EZ\{0} (3.19) 

corresponding to charged fields of charges {qi} inserted at points {xi} and local 
observables located at {Cj}, with Xi and Cj in the positive time lattice. 
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From (generalizations of) the lower bounds (3.16) it follows that ||Te|| 7^ 0, i.e. 

the states (3.19) have finite energy. 

One can define non-local charged fields by 

$(x , q, E) \xi, gi, ...x n , q n , Ci, C m > E = \x, q, xi, qi, x n , q n , Ci, C m >£ 

(3.20) 

with x° < x\ < ... < x°. 

We now consider the region of coupling constant space where all correlation func- 
tions of non-zero total charge vanish. Then H(E) decomposes into orthogonal 
sectors labelled by the total electric charge q: 

H{E) = ® q H q {E). 

Cluster properties show that the sectors 7i q (E), q 7^ 0, are charged sectors, in the 
sense described in sect.3.1. 

From the above construction it is clear that, a priori, the charged sectors 7i q (E) 
depend on the choice of the distribution E. It is natural to ask if 7i q (E) is 
orthogonal to H q (E'), for E 7^ E' . In order to answer this question , one considers 
the generalized correlation functions (3.16). 

One can define a scalar product between states in H q {E) and H q (E') by 

E <x, q\x', q > E '= G 2 ,o(rx, -q, E, x', q, E') 

From generalizations of the upper bounds (3.19), (3.20 ) it follows that H q (E) _L 
n q {E') if 

(^(E — E'), A~ 1 (i? — E')^j (3.21) 

diverges and one easily realizes that this divergence occurs if E and E' do not have 
the same "behaviour at infinity". 

There is an interesting choice of E and E' which naturally leads to a divergence 
in (3.21). For example, if one chooses E to be supported in a spatial cone S 
with apex in and opening solid angle less than 4n, the states obtained via O.S. 
reconstruction are the lattice approximation of the states discussed by Buchholz 
[26] in the algebraic approach to Q.E.D.. If the field E' is chosen to be localized 
in a disjoint spatial cone S 1 ', then (3.21) diverges. In particular, if we choose S' 
to be the cone obtained by a rotation of S then this divergence shows that in the 
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continuum limit (if it exists) the rotations cannot be unitarily implemented on 
"Buchholz states". For more details see [2,3,19]. 

We conclude with some remarks about particle structure analysis on charge ±1 
sectors. By ispection of the proof of Theorem 3.2 (see [2,3]), one can argue that 
for \x — y\ large : 

G* (x, -1, y, 1) ~ <P*v) J ^en(i)e ,(A '"- + ^ (3.22) 

Pxy '-^Pxy — &x &y 

where p xy is a current of flux 1 and connected support E = E(x) — E(y), z(p xy ) 
is a statistical weight and dp ren (A) is a positive measure of the form 

dfiren(A) = -^—d^ ren {A)e x ^ (3.23) 

In (3.23), d/j,p reri (A) is a gaussian measure with mean and covariance (3 ren (5d)~ 1 
(+ gauge fixing), where (3 ren is a renormalized coupling constant [(3 ren {K,f3) = 
(3 + 0(c(k))] and X{dA) is a sum of gauge-invariant "irrelevant" terms, in the 
jargon of the renormalization group. Hence the large distance behaviour of Gfjcb 
which is independent of the irrelevant terms, X, should essentially be given by a 
product of two factors: one is due to the fluctualing current line p xy , and, from 
the analysis in terms of excitation expansions of [23], it is expected to produce an 
Ornstein-Zernike decay 

| a; _ y |-3/2 e -m( K ,/3)|x-y| 

corresponding to a particle of mass m(«, (3) ~ — lnc(«) +0(f3). The second factor 
can be argued to contribute another power correction to the exponential law 

exp [-^(£, A" 1 ^)] ~ \x - y\- c ^\c > 

It is due to the soft photons accompanying an infra-particle. Therefore, as x° / 
oo, one expects that 



< 0(0, ±1, E)Q, (f)(x, ±1, E)Q > 

x° / oo (3.24) 



( x 0^3/2+cf3 ren 

Equation (3.24) exhibits the infraparticle nature of the charged particles in the 
Higgs model. In fact, the vacuum expectation value of a charged field vanishes 
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in the region of coupling constant space considered, and a comparison with the 
general formula (3.5) shows that the mass operator M, as defined in (3.4), does 
not have a discrete eigenvalue corresponding to a sharp one-particle state. 

Remark 3.3 (Q.E.D.) In lattice Q.E.D., the electron field is defined by 

and the conjugate euclidean field is 

*(x,-l,E) = q x -f e- t{A > E(x)) 

Typical local observables are Wilson loops and string fields. Correlation functions 
of charged fields and local observables satisfy the hypotheses of the Reconstruction 
Theorem 3.1 and, for k sufficiently small, bounds analogous to those in Theorem 
ce[5 3.3 (as discussed in [2]). 

This permits one to reconstruct non-local electron-positron field operators 
ty a (x, q, E), q = ±1, which are defined by 

V a (x,q, E)\xi,q!,ai, x n q n a n , d, ...C m > E = 
= \x,q,a, x^qxax, ...,x n ,q n a n ,C\, ...,C m > E 
for x = ( 2 , x) , x 7^ Xj . 

The field ^ a (x,q,E) can be viewed as the lattice approximation of the formal 
operator (1.1) introduced by Dirac. An analogue of equation (3.24) is expected to 
hold for k small enough, on the basis of the proof of bounds on electron-positron 
correlations. It would exhibit the infraparticle nature of electrons and positrons. 

3.3 Monopoles and duality 

Model B, in Villain form, at k = is dual to model A in the limit k S oo 
and it describes an interacting theory of "photons" and Dirac monopoles. Since 
monopoles can be viewed as solitons in such a model, we may appeal to the strategy 
outlined in sect. 3.1 to construct a monopole field operator by introducing disorder 
fields, whose expectation values are the Euclidean Green functions of the Dirac 
monopoles [2,3,19]. For this purpose we introduce a real 3-form B on the lattice 
Z 4 , given by B = 2ir*E and define the disorder field by 

^(x 1 , ?1 ,...,x n , ?n ,S) = e-^ll^+ 5A - 1 ^--)ll 2 -ll^ll 2 > (3.25) 
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where x { G Z 4 , q t G Z\{0}, 



i 

and a; is a 27rZ-valued 3-form satisfying 

d(B -oj) = 0. (3.26) 

We now explain why expectation values of D u are correlation functions of monopoles. 
First, we use the Hodge decomposition to rewrite 

n = dA- 1 5n + 5A- 1 m (3.27) 

where 

m = dn (3.28) 

Let n[m] be an integer-valued solution of the cohomological equation (3.28). Then 
every other solution is of the form n[m] + d£, where I is a to 27rZ- valued 1-form. 
We define a real- valued gauge field A by setting 

A = e + £ + SA- 1 n[m] (3.29) 
and one can easily verify that 

(D UJ (x 1 ,q 1 , ...,x n ,q n ,B)) = 

l^m:dm=0 J 1 l<ij> u/± <tJ > e , , 

v rrr ha P -§\\dA+sA~^m\\ 2 ' ' ' 

2^m:dm=0 J 1 l<ij> aA <iJ> e 2 " 

A is the euclidean "photon field" , and the Hodge-dual of (m — uj) is supported 
on a set of lines in the dual lattice Z\, which can be interpreted as Euclidean 

2 

worldlines of Dirac monopoles; uj itself can be viewed as a Dirac string if we take 



= Uxi + Uc 



where *uJ Xi has support in an open line at constant time with one end at X{ and 
the other end joining a compensating current uJoo at infinity (monopole b.c). 
Then B(xj) — uj Xi is exactly the lattice approximations of the magnetic field of a 
monopole located at Xi, together with its Dirac string uj Xi . 
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In the presence of the disorder field the Euclidean observables of the model 
must be modified so that expectation values do not depend on the choice of the 
"Dirac string" to. In particular the Wilson loop 0(C) is now replaced by 

O u (S) = 0(C) Jl e~ l{SA ~ luj) (3.31) 
P es-.dS=c 

where S is a surface; see [19,2]. 

Correlation functions in the U(l) gauge theory are then defined by 

m 

; Si, S m ) = (D u yXi , q\ , x n , q ni 

B)l[O u (Si)) (3.32) 

i=l 

By duality e.g. 

n 

(D u (x 1 ,q 1 ,...,x n ,q n ,B))= lim (TT Qfa, ft, E)) K (3.33) 

1 = 1 

where (-) K denotes the expectation value in the Higgs model A and E =* B. 
O.S. reconstruction theorem applied to G^ m provide us with non-local Dirac 
monopole fields Ai(x, q, B) and, for (3 large enough, Dirac monopole sectors 7i q (B). 
Furthermore a particle-structure analysis along the line of sect. 3.2 exhibits the 
infraparticle nature of the Dirac monopole of charge ±1. 

Remark 3.3 The monopole construction outlined above can be applied to every 
{/(l)-gauge theory without matter fields. This constructions and variants thereof 
have been applied in numerical simulations of lattice theories in [27], with the aim 
of detecting phase transitions between phases corresponding to different behaviour 
at large distances of the monopole Green functions, which has been described in 
the dual picture in Theorem 3.3 and Remark 3.2. 

4. Charges and monopoles 

We start this section with an outline of the new problems appearing in an at- 
tempt to extend the construction sketched in the previous sections to models 
where charges and monopoles coexist (class B,C). 

4.1 Dirac quantization condition 

We have seen that, for k = 0, model B (in Villain form) can be written explicitly 
in terms of a "photon field" A and a "monopole field" *m. The Dirac strings of 
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the virtual monopole loops in the partition function sweep out surfaces described 
by the support of *n. The term Si in the action introduced in eq.(2.6) couples the 
gauge field 9 to a charged matter field. In our expression for the partition function, 
the charged matter gives rise to charged loops describing the worldlines of virtual 
particle-antiparticle pairs. We must ask whether the statistical weights of these 
charged loops depend on the location of the Dirac strings of virtual monopoles. 
To answer this question, let us rewrite the action (2.6) in terms of A,m,(p,£. We 
obtain, for q = 1, 

S = -\\dA\ | 2 + |(m, A _1 m) + -\\d<p + A - SA-^m] + t\ | 2 (4.1) 
2 2 2 

Independence of the Dirac strings corresponds, as recalled above, to independence 

of the choice of the two-form n[m] satisfying dn[m] = m, as in (3.28). 

Let ( be a Z-valued 1-form. Then by Poisson summation formula, 

^ e -§ Wdip+A-SA^nlml+eW 2 = ^ e ~ sfc I Kl \ 2 ^(Cdf+A-SA^nlm]) 

e C 
Integrating over <p we obtain 

5( = 0. 

We conclude, using the Poincare lemma, that there exists a Z-valued 2-form, x, 
such that 

C = s x , 

so that 



(C,5A 1 n[m]) = (Sx, 5 A 1 n[m]) = (x, dSA 1 n[m]) = (%, n[m]) — (x, 5 A 1 m) 

where the last equality follows from Hodge decomposition. 
The second term depends only on m and 

e »fc«W) = i 

This is nothing but the Dirac quantization condition, since the electric charges 
appearing in the partition functions are q e e Z and the magnetic charges are 
q m £ 27rZ, so that 
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q e q m e 2ttZ (4.2) 

From the above proof it is clear that if we consider the expectation value of e l ( d,E \ 
where the "electric current" 1-form E is not integer valued, as introduced in sect 
3.2, we would encounter an inconsistent dependence on the choice of Dirac strings. 
This problem is often neglected in the physics literature, wher virtual monopole 
loops are sometimes ignored, assuming that monopoles are "very heavy" , see e.g. 
[28]. But, in order to arrive at a fully consistent definition of gauge-invariant 
charged field Green functions, one must cope with it. 

The natural suggestion is to replace the electric field E(x) of sect. 3.2 by an 
integer-valued electric current, a "Mandelstam string", starting at x and ending 
at the location of some compensating charge which will eventually be sent to 
infinity. 

A naive idea would be to simply replace the variable <j) x e % ( A,E ( x '\ used in the 
construction of charged states in model A, by 

c R e i<P* e i(e,7?) (4.3) 

in model B, and then take the limit R S oo. In (4.3), 7^ is a unit 1-form with 
support on a straight line in a fixed-time plane from x to some point at a distance 
R and then joning that point to a fixed point in the time-zero plane, and c x is 
some normalization factor. 

Consider the 2-point function: Integrating out ip we obtain 



( e i<P* e M>r?) e -i<Pv e -M>r?)) = 2_ J JJ ^ <lJ>e - s «W ^ z{p)e^ e ^-<^ 

<ij> p:Sp=8 x -S y 

(4.4) 

where p are Z- valued currents whose statistical wheight is denoted by z(p), and 
Z is the partition function. The phase factors appearing in (4.4) define Wilson 
loops, and, in the region of coupling constant space where monopoles are particle 
excitations, the expectation value of the Wilson loop is known to exhibit perimeter 
decay. Hence we expect (4.4) to vanish exponentially fast as R S 00. 
This decay is dominantly due to the self-energy of the strings 7. This effect could 
eliminated by adjusting c. But we would then be left with an interaction term 
between the strings which, being attractive and extending over the full string, 
tends to infinity in the limit R S 00, and this appears to render the renormalized 
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Green function divergent. Furthermore, since the interaction term depends on the 
time distance between strings, it cannot be renormalized away without violating 
O.S. positivity. This problem could be solved if we replace a "straight Mandelstam 
string" by a sum over "fluctuating Mandelstam strings", weighted by a measure 
which is concentrated on strings fluctuating so strongly that, with probability one, 
their "interaction energy" remains finite in the limit R /* oo. 
A proposal for a "natural" measure can be inferred from a representation of cor- 
relation functions of a U(l) scalar field, x, with coupling constant (3 X , coupled to 
an external U(l) gauge field 6>, in terms of random walks. E.g., in d dimensions, 
with <•>(#) denoting the corresponding expectation value, 

<e*-e-'*>(*) = E {2d) ~^ vhl Z ^jM (4.5) 

where u xy is a path ("string") from x to y, \uj xy \ its length, Z{9) is the partition 
function of the system, and Z(9\u xy ) is a partition function modified by u xy , see 
[29]. 

In d > 3 dimensions, for j3 x sufficiently large and 9 = 0, the string u xy is known 
to be rough. In fact 



{e lXx e lXv )(0) — > const 

as \x — y\ — > oo, and it is believed that it remains in the rough phase for a 
class of positive measure of L r (l)-external gauge fields, (roughly speaking, if d6 is 
sufficiently small in average). 

This suggests that two-point functions of a U(l) scalar field might yield an ap- 
propriate measure on "Mandelstam strings" , if j3 x is large enough. 
For a better understanding of our construction in the compact models, we show 
how to reproduce the results already obtained in the non-compact models following 
the above ideas. 

4.2 The non-compact model revisited 

Let T(R,t),R G Z + ,t G Zx/2) be a cube of height 2\t\ in the time (x -) direction 
and with sides of length R in coordinate planes x° =const and centered at the 
origin; see fig. 2. We define A(R,t) to be given by dT(R,t) fl {x : x° ^ 0}, for 
t ^ 0.(To simplify our notation, the restriction of a lattice operator to A(i?, t) will 
be denoted with a subscript A, the specific sublattice which we are referring to 
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being identified from context.) We introduce a real-valued 0-form A on A(R,t), 
with action given by 



S A (X) = ^\\d A X + A\\l (4.6) 

where (, )a is the inner product in A(i?, t), and ||/|| A = (/, /) A - 
Denote by (■) a(r,x°)(A) the (normalized) A-expectation corresponding to the ac- 
tion S A (X). Set R± = (±±,R, 0, 0) and, for x° > 0, y° < 0, define 



G(x,y) = Jim ((e^e-^ + ) A(i?)X o ) (A)(e-^e a «-) A(i? , y o ) (A)e^<« + «->$ x $,) 

(4.7) 



where (•) is the expectation value of model A. 
By explicit computation 

( e iXx e~ iXR +) A (R x°){A) = e ~2WI ( - ( - Sx ~ 5R +^ A A 1 ( S ^- 5 R+y) e i(E A (x)-E A (R + ),A) ^g-j 

where 

E A (x) = d A A^8 x (4.9) 

Hence 



G(x, y) = lim ([e l ^ A( ^- EA( ^' A) $ x $J[e l(EA(i? - ) - SA(i? + )+5 <«+«->' A) ]) 

R^oo L J 

= c ( e i(E(x)-E(y),A)Q x Q^ ( A1Q ^ 

where 5<,r +j r_> denotes the 1-form with support on the link < R + R_ >, whose 
value is 1 on that link, and 

c= lim ^ e i ( £; A(i?-)-i?A(i?+)+5<H + R_>,A)^_ 

R^oo 

To prove (4.10), we note that the form E A (x) — E A (y) decays like d~ 3 in the limit 
R /" oo, and E A (x) — > E(x), where E is given in (3.8). 

We observe that the two-point function of the auxiliary matter field e lX exactly 
reproduces the exponential e l ( A,E ^ needed in the construction of charged states. 
Furthermore our construction respects O.S. positivity and, when the limit R S oo 
is taken, lattice translation invariance is restored. 
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Sectors corresponding to different choices of the one-form E can be reproduced 
by modifying S\ (A). 

4.3 Charged Gelds in compact models 

It is natural to generalize the construction of the previous section to the compact 
model B. 

We introduce a {7 (1)- valued scalar field x and a 27rZ-valued 1-form r on the lattices 
A(i?, t) defined in the previous section. An action functional for x and r is defined 
by 

SA(x,r) = ^\\d AX -e + r\\i (4.12) 

We denote by (-)A(R,t)(9) the expectation value w.r.t. the action (4.12). 

With notations analogous to those of the previous section, we propose the following 

definition of the two-point function of the charge-1 field in model B: 

G(x,y) = 

jgm c£ic*i ( (e-*- ) A(fliX0) (9) ( e *. e -*«- ) A(fl>y0) (0)e ie <«+«- > e*- ) 

(4.13) 

(See fig. 2) In (4.13) c% is a normalization factor given by 

<£, q = i ( (e~ iX * + d XR + > A(H,i/2) (*) e-««- ) A(fl ,_ 1/2) (^)e"^<« + «- > e*<- + -- > ) }"* 

(4.14) 

where q is an integer and x± = (±|,x). An approximate evaluation of (4.14) 
suggests that c^ q might be bounded in R (see Appendix A). 
The charge-1 two-point function (e lXx e~ lXn + )a(r,x°)(9) is periodic in 9, with pe- 
riod 1. Hence, it has a Fourier representation 

(e^e-^+^o)^) = J d m+ { 3x )e^^ (4.15) 
where j x are Z-valued one-forms satisfying 

Sj x = S X - S R+ 

and dfiR + (j x ) is a complex measure on the space of 1-forms j x . One can integrate 
out the matter field <p in (4.13) and express the contribution as a weighted sum of 
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Z-valued 1-currents. In the numerator of (4.13), a contribution due to a 1-current 
u satisfying 



5u = 5 X — 5 y 

will appear. As a result of this representation of expectation values in terms of 
sums over integer-valued currents, equation (4.13) yields a representation of the 
Green function G(x, y) as an expectation value of a weighted sum of Wilson loops 

W u (j x ,j v ) = e t{j * +J y +s < R + R -> +UJ > d) (4.16) 

(See fig.2) 

To analyze (e~ lXx e lXR + )a(r,x°) (0)> we perform a change of variables analogous to 
(3.36). Define 

v = d A r (4.17) 

and let r(v) be a fixed 2nZ- valued solution of equation (4.17). We then define a 
real-valued scalar field A by 

X = X + S A Al 1 r[v] 
Changing variables from x, r to A, v we obtain 

(e-^e^ + > A( ^ O) (0) = 

E„-^=o/n l ^ e "^ IMAA+5AA ^ lu+e|l ^ ei(A ~ <5AAXlrM ' ,5:c "' 5i?+)A 

= e -^((^-^ + ),A^ 1 (^-5 K+ )) e _ j(i , A(:l;ii?+) , 9)F(£;A(x ^ +)| ^ ) (4 _ lg) 

where 

V p -^{v+d6,A^{v+d6)) A -i(E A ,r[v]) 

F(E\d9) = ^"=^=o c _ (4.19) 

V p -^(v+dd,A]i 1 (v+d0y)A 
l^v:dv=O fi 

and 



E A (x,R+) = E A (x) -E A (R+). 
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The action (4.12) describes a U(l) spin system in the presence of the external 
gauge field 6>; the first two factors on the r.h.s. of (4.18) correspond to the spin- 
wave (gaussian) approximation of the correlation function, the last one to the 
contribution of the vortices, described by v, to the spin system. 
For (3 X sufficiently large and dO = 0, the spin system is known to be in a phase 
with long range order, where the spin-wave approximation is known to capture 
the essential physics at large distances. It is argued in [29], as remarked above, 
that, in the presence of the gauge field, the system remains in a phase with LRO, 
provided dO is sufficiently "small" in average. Hence, for f3 x ,f3 sufficiently large, we 
expect the vortices to form a dilute gas, so that the factor F in (4.18) is expected 
to represent a "small" correction to the spin-wave approximation. However, the 
presence of F is crucial for the periodicity in 9 which ensures independence of the 
Dirac strings of virtual monopole loops. 

Following the arguments of sect. 3.2 (but see also [2,14]), one can argue that, for 
k small enough, the contribution of the matter field at large distances essentially 
yields a renormalized U(l) gauge theory with a coupling constant 

(3' ren =P + 0(K 4 ) 

which, according to [12,24,3], and sect. 3.4, is expected to renormalize to a non- 
compact gauge theory with coupling constant 

/3ren = Pren ~(~ e ^ren) 

in the scaling limit, provided (3' ren is large enough. Accordingly, the 2-point func- 
tion is expected to behave like 

G ( Xf y)„L J VAe -MdA\\* ^ z{pxy y(A, Pxv+ E( X )-E(y)) 

Pxy Pxy — &x &y 

F G (E(x)\dA)F G (-E(y)\dA), (4.20) 

at large distances , where F G is a gaussian approximation of F estimated in 
Appendix A by 

F G (E\dA) ^ e -le-° (Px) \\E\\ 2 +e- 0( ^ ) \\E-5 A Al 1 dA\\ 2 ] ^ ^ 

Appealing to universality in the scaling limit, we believe that the large distance 
properties of F are indeed correctly captured by the gaussian approximation, F G . 
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Green functions with vanishing total charge G m+n (xi, qi, x m , q m , y\, q[, y n , q' n ), 
> 0, yf < 0, are defined by 

G m -\-n(xi, qi, x m , q m , yi,qi, y n , q n ) 
■ 

3 = 1 

n 

II <U ( e ~ WkXyk eWkXR - )A(R, y o ) (0)e iq '^ e^<«+«- > ) , (4.22) 
fc=i 

where q = Xw=i 9j an d c <? 9 i s as i n (4-14). 

The large distance behaviour described by (4.20) implies cluster properties en- 
abling us to reconstruct charged sectors and non-local charged field operators 
${x,q,E) for model B. Equation (4.20) also suggests that the particles of charge 
±1 are infraparticles (see sect. 3. 2), since the contribution due to F G does not 
change qualitatively the large distance behaviour [30] . The same ideas also apply 
straightforwardly to model B with q ^ 1. 

These observations suggest that the Green functions of charged fields in the present 
model have large-scale properties similar to those in gauge theories without dy- 
namical monopoles, constructed according to Dirac's original proposal, provided 
(3 X is chosen large enough. On short distance scales, the Dirac quantization con- 
dition implies integral quantization of electric flux accompanying a charged field 
<&(x, q, E) in all models with dynamical monopoles, while electric flux is not quan- 
tized in Dirac's construction. But, on large distance scales, the constraint of in- 
teger electric flux becomes "irrelevant", (i.e. this constraint renormalize to zero). 
The phenomenon described here is analogous to the phenomenon of symmetry 
enhancement studied in [12,31]. 

Remark4.1 The construction described above can be generalized to Green functions 
of charged fields in deconfined (non-abelian Coulomb) phases of non-abelian gauge 
theories with some gauge group G, by simply replacing the auxiliary two-point 
correlation functions of the U(l) scalar field x coupled to the U(l) gauge field by 
two-point correlation functions of a non linear sigma model with target G coupled 
to the gauge field. These correlation functions can be expressed [29,32] in terms 
of sums over path-ordered exponentials of the gauge field. It is useful to recall 
that an inequality analogous to the one quoted after (4.5) holds. If G contains a 
U(l) subgroup and the field, g, of the non-linear sigma model, transforms under 
a representation, U, corresponding to a character non trivial on this subgroup, 
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then, at weak coupling, the two-point function in an external gauge field of zero 
curvature, (L r ((7 x )L r ((7~ 1 ))(0) is bounded uniformly from below in \x — y\. One 
expects that this long range order persists in the model coupled to the G-gauge 
field, provided the curvature of the gauge field is small enough in average [29]. 

4.4 Monopoles in compact models with matter 

In this section, we propose a construction of monopole Green functions in compact 
models with matter (class B) obtained through a duality transformation. 
First we identify the model dual to B. It can be described in terms of a Z- valued 
2-form 7 and a Z-valued 1-form a (defined on the dual lattice) with an action 
given by 

S(an) = h\da - qi\\ 2 + ^\\d^\\ 2 (4.23) 
We can replace a by a real-valued 1-form A by inserting the constraint 

m p 

where m is a Z- and p is a 27rZ-valued 1-form, and we have applied the Poisson 
summation formula. 

The partition function of the dual model can then be written as 

Z = J2 f J! dA <l3> e-^ dA -^-^ d ^ el(A ' P) ' ( 4 - 24 ) 

[7] <ij> p:Sp=0 

where we exploit the gauge-invariance of the A measure to impose the constraint 
dp = 0. 

The interpretation of formula (4.24) is clear: it describes an abelian gauge theory 
with charged matter fields and dynamical monopoles satisfying Dirac's quantiza- 
tion condition. The 1-form p describes the worldlines of electric point charges (with 
values in 2-7rZ), the 1-form *d-) describes the worldlines of magnetic point charges 
(which are integer-valued). This interpretation is consistent with the feature of 
duality transformations that they exchange electric and magnetic charge. 
Accordingly, a Green function of charged fields in the dual model corresponds to a 
monopole Green function in the original model. If we try to follow the construction 
of the previous section we encounter the following problem: we can reabsorb the 
"Mandelstam strings" appearing in an expansion of the x~correlation functions, as 
in (4.5) and (4.13), in a shift of the current p. One may say that these "Mandelstam 
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strings" are "screened" by the charge loops appearing in the dual model. Starting 
from the action (4.23), we meet the same phenomenon in the following way: Since 
the gauge field a is integer-valued and the electric charge of the auxiliary field x 
is 27r, the dynamics of x is independent of a, (i.e., coupling x to a does not have 
any physical effects). 

A way out of this difficulty appears when q ^ 1, because the magnetic charges 
appearing in the model are multiples of q, so that one can introduce electric loops 
of charge ^l, still satisfying Dirac's quantization condition. 

Hence, instead of inserting correlation functions described as sums over a single 
fluctuating string of charge 27r, one can introduce correlation functions described 
as sums over q distinct fluctuating strings of charge ^ which are not screened. 
For this purpose we introduce a scalar field x with value in [—nq, irq] and a 2tvZ- 
valued 1-form r with action 

S,(x,r) = ^\\ d -^^ + r\\l (4.25) 
z q 

and denote the expectation value w.r.t. the action (4.25) by (•)A(i?,t)(A), (see sect. 
4.2). The charge-f two-point correlation function (e %Xx e~ %XR + )a(r,x°)(A) has a 
Fourier representation 

(e ix -e- ix ^> A (fl,«o)(i4) = J d m+ ( 3x y^ (4.26) 
where j x are -valued f -forms satisfying 



5j x = 2n5 x - 2n5 R+ 

and d(iR + (j x ) is a complex measure on j x . As a 2-point function of charged fields 
in the model dual to B, for q ^ 1 , we propose to consider 



G*(x,y) = ^cfcf((e-^e-^ + ) A(i? , :c o ) (a)(e^e-^-) A(i? , J/ o ) (a))* = 
fl ln^c^((e^e-^^ 

(4.27) 

where (•)* is the expectation value of the dual model, and w is a 27rZ-valued 
1-form satisfying 



Suj = 2ti{5 x — S y ) 
30 



(4.28) 



(See fig. 3) Note that G* is independent of the choice of u satisfying (4.28). 

The corresponding two-point function in the original compact model is obtained 

by applying an inverse duality transformation, by replacing the Wilson loops 



W u (j x ,j y ) = e ^+Jv+^ <R+R _ >+ u,A) 

appearing in an expansion of (4.27) (as discussed in the previous section, see (4.16)) 
by 't Hooft disorder loops. 

To construct the 't Hooft loop, we have to identify a surface whose boundary is 
given by the loop L = j x + j y + 2n5 <R+R _ > +* uj. Let us denote its ^^-valued 
Poincare dual by H(j x ,j y ;oj). The 't Hooft disorder field, Du(j x , jy), is obtained 
by shifting d0 by E in the original model, i.e., 

D u (j x ,j y ) = e -^(W M +n^,3y^)+n\\ 2 -\\de+n\\ 2 ) 

Combining (4.26) and (4.27), the two-point monopole function in model B is given 
by 

G*(x,y) = lim c^c^i I dfi R+ (j x )d m Jj y )DUj x Jy)) (4.29) 

R— >oo J 

Equation (4.29) expresses G*(x,y) as an expectation value of a weighted sum of 
't Hooft disorder fields D u (j x ,j y ). In order to better understand the meaning of 
(4.29), notice that, as in sect. 4. 2, one can easily show that 



(e i x*e- i *«+) K{R , x o ) {A) = 



e -^((^-5H + )-AI 1 (^-^ + )) el ( EA (x,i? + ),A) F(EA(x ^ +) |^ rfA) _ (4 _ 30) 

From (4.26) and (4.30) it follows that dji{j x ) behaves, for large f3 x , as an approxi- 
mate Dirac 5 measure around E\(x, R + ). This, in turn, implies that dE(j x , j y ;u>) 
is peaked around 

B A (x,R+) +B A (R-,y) +* oo + 2n*5 <R+R _ > , with B A =* E A 

Here we see the dual version of the phenomenon occurring for charged fields: the 
magnetic current appearing in the definition of monopole Green functions is forced 
to be 2-7rZ valued so as to respect Dirac quantization condition. But its average, 
at sufficiently large scales, approximates the ordinary magnetic field £?, of the 
classical monopole. 
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The deviation from a 5-measure of dpR + (j x ) is "small", for j3 x large. This may 
explain why, in numerical simulations [33], a naive application of the definition 
of monopole Green functions in theories without matter fields appears to give 
reasonable results. 

Redefining j3 x = q 2 /3 x , keeping (3 X fixed and taking the limits g / oo, k \ one 
recovers the monopole Green functions discussed in sect. 3.3. 
Using a Peierls type argument and renormalization group techniques discussed in 
[2,3,12,24], we may argue that the monopole Green functions in a theory with 
matter fields have the same qualitative large distance behaviour as the monopole 
Green functions in a pure U(l) gauge theory, if f3 and k~ 1 or q 2 / (3 are sufficiently 
large. (In the dual model the renormalization of the electric currents p requires j3 
to be large, and the renormalization of the monopole currents d^y requires k~ x + 
0(q 2 //3) to be large.) 

4.5 Dyons 

Model B in Villain form is self-dual in the limit k / oo, where it becomes a Z g 
gauge theory. If a topological term (2.9) is added one can show, see [11], that 
a modified duality transformation is an approximate symmetry. When combined 
with the 27r-shift of © it generates an SX(2, Z) duality group of approximate 
(large-scale) symmetries of the model. 

Here we briefly describe how to derive the duality symmetry and then discuss its 
action on charged and monopole Green functions. Green functions of dyon fields 
are obtained in a natural way. 

A wedge product on the lattice compatible with O.S. positivity can be defined [17] 
as follows: Label a (positively oriented) /c-cell in the lattice by a site, x, and a set 
of k directions fi = {//i, ...pk}', i-e., 

k 

c k = (x; fi u /i k ) = x + i^ e ^i > ( 4 - 31 ) 

i=i 

where £ Mi G [0, 1] and e Mi is the unit vector in the \ii direction. 
Given a p-form A and a (/-form B we define 



AAB(x;p) = Y ^' ->fh) B ( x + e m + - + e Mp ;/V+i,...,^ P +q) 

+ ]-A(x] rn u r/i p )B(x + e rfll + ... + e rfip ] r/i p+1 , r(i p+q )} (4.32) 
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where r inverts the time direction, leaving the other ones unchanged. 
It is easy to verify that 



d(A A B) = dA A B + (—1) P A A dB (4.33) 

but in general 

AAB ^ —(—) p+q B A A 
Furthermore, on the lattice, we do not have an analogue of the continuum identity 

J AAB = (A,* B) = (*A,B) (4.34) 

for p + q = 4, where * is the Hodge dual and (,) denotes the inner product on 
forms in R 4 . 

From the definition (4.32) it follows, however, that we can define two modified star 
operations, A \— >* A, A \— > A*, such that for p + q = 4, 

J2(A A B)(c 4 ) = (A* B) = (A*,B) (4.35) 

c 4 

and with the property that if 6* A = 0, then 5* A = = 5 A*, but A* ^ A and 
A** ^ A** A^ A. Furthemore, if A is Z-valued, A* and *A are Z/2-valued. 
Using the rule (4.33) one can rewrite the topological term (2.9) in terms of the 
magnetic currents m = dn : up to boundary terms vanishing if O-Dirichlet b.c. 
are imposed on dd and n, it can be rewritten as 

i-^q^2(mA6 + 6 Am + nAn)(c 4 ) (4.36) 

c 4 

Its meaning is clear: the first two terms associate to every magnetic current m a 
line of total electric flux J|g along the support of m* and *m (Witten effect [34]); 
the last term is the self-intersection form of the surface *n. 
To discuss the duality group it is convenient to define 

A±A* 

for each /c-form A, and to introduce, following [11], the complex coupling constant 
parameter 
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(and its complex conjugate, £) . 

The pure gauge action of the model can be rewritten as 



So(e,n) = ^(\\(de + n) + \\ 2 - licil(^ + n)_|| 2 (4.38) 

Z7T Z7T 

and, integrating over the matter field, the partition function is given by 

z = J2 I II de <lJ>e - s ^ eMe,p) ( 4 - 39 ) 

n J <ij> p:Sp=0 

The SL(2, Z) group acts on £ by 

Its generators are called S : C — > and T : £ — > £ + 1. 
We now show why S is an approximate symmetry of the partition function. 
Introducing real 2-forms A, A, denoting by [a] the equivalence class of Z-valued 
two-forms, i.e., [a] = {a' — a = 5V, V(cs) G Z}, and performing the change of 
variables {6>,n} — > {A, m}, as in previous section, one obtains that 

z= Yl f U dA <^>U dX pU d ' x p 

m:dm=0 [a] <ij> P P 

e -^\\M\\+^\W\\ l (dA+n[rn], ( 4±+ ( 4 r+q *) 

Integrating out A we obtain the constraint 

If the identity A** = A were true then it would follow from (4.40) that there exists 
a real- valued 3-form, £, such that 

A =(<?£- 27rg<7)+ , A = (^ - 2nqa)- 

Formally, A* and A* have the same continuum limit. On large scales, we may 
therefore replace * by *. Then A + and A- are just the selfdual and anti-selfdual 
components of A and, indeed, A** = A. 
So, within this approximation, 



E E 

m:dm=0 [a] 



J c 

(4.41) 
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Passing to the dual lattice and introducing a [/(l)-valued 1-form 9 and a Z-valued 
1-form I and setting 



i = q*6 + 2nq*£, 

we can rewrite (4.41) as 

Z ~ E E/ II de<y>e-^ll^- 2w '^+ll^ + ^ll^- 2w, ^-ll^^^^ 

m:dm=0 o <ij> 

(4.42) 

Equation (4.42) proves that the £ generator SX(2, Z) induces an approximate 
symmetry exchanging 

C -> - m -> -2tt p p -> — . 

C ^ 

Analogously, the T-generator induces an approximate symmetry exchanging 

as one can easily see from (4.15) (4.36); the symmetry becomes exact for T 2 since 
^ and are actually Z-valued, as the original p current. 

Our construction of the charged 2-point function, G(x,y), is based on the intro- 
duction of a weighted sum of Wilson loops 

W u (jxJ y ) = e iqie '^ +J y + * UJ+5 < R + R '> ) (4.43) 

where a; is a current line from x to y, and our construction of the monopole 2-point 
function G*(x, y) is based on the introduction of a weighted sum of 't Hooft loops 

yjx i jy ) 

e - l £{\\(d0+n+XU x ,j y ;co)) + \\ 2 -\\(de-n) + \\ 2 } e i£{\\(d0+^^ 

(AAA) 

with the notations of previous sections. 

The above discussion makes it clear that the S generator of SL(2, Z) approximately 
exchanges the two correlation functions in the dual models, while the T generator 
acts (approximately) by multiplying the 't Hooft loop (4.44) by the Wilson loop 
(4.43), raised to the power —2, hence producing the 2-point function of a dyon, 
whose electric charge is — 2q at = 0. 
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5. 't Hooft-Polyakov monopoles 

5.1 Preliminaries 

In order to explain our proposal for the construction of Green functions for the 
quantum 't Hooft-Polyakov monopole, it is useful to recall the definition of topo- 
logical invariants characterizing the classical 't Hooft-Polyakov monopole in the 
Georgi-Glashow model and then to exhibit their lattice counterparts. 
The field content of the classical Georgi-Glashow model consists of a connection 
A on a principal 5 , L r (2)-bundle, Psu(2), on R 4 and a scalar field $ with values 
in su{2) = Lie{SU{2)), which is a section of a vector bundle associated to Psu{2) 
with fiber su(2) equipped with the adjoint action of SU{2). 
The Lagrangian density of the model is given by 

^=^l^| 2 + ^|^| 2 -^(|$| 2 -l) 2 (5.1) 
where Fa is the curvature of A, and Da is the covariant derivative 



D A $ = d$ + [A, $] 

Moreover | • | denotes the Killing norm on su(2). Let us assume that A and $ verify 
appropriate regularity conditions, and |1 — |<&(a?)||, \Da$>\(x), Fa(x) have decay 
properties at oo discussed in detail in [35]. Then, at fixed time, a finite-energy 
static configuration (A, $), with A = 0, defines a homotopy class, [(A, <&)] G 
TT2(SU(2)/U(1)) ~ Z, labelled by an integer topological charge 

N = -i- / Tr(F4 A D A §) = ±- [ Tr($F A ) (5.2) 
An J R s An J ^ 

where is the 2-sphere at infinity. Furthermore, the field $ defines a homotopy 
class [<&] G n 2 (S 2 ) which can be characterized by the Kronecker index 

N = —[ Tr($Arf$Ad$); (5.3) 

with $ = $/|$|, and Sr = {\x\ = R} provided R is sufficiently large. Then (5.3) 
is independent of R, and [<fr] = [(A, <&)]. A solution of the equations of motion 
with 7^ is called a classical 't Hooft-Polyakov monopole. 
The magnetic charge of a Dirac monopole in R 3 is given by the first Chern num- 
ber of a ?7(l)-bundle over a topological 2-sphere containing the monopole. [This 
geometrical definition differs from the one used previously by a factor l/(2n).] We 
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show how one can associate to a configuration (A, <£>) and a cube c in R 3 a £7(1)- 
bundle in such a way that the first Chern number of the bundle can be identified 
with the magnetic charge contained in c. We then relate the magnetic charge to 
the topological charge previously defined. 

We choose a cube c in R 3 ; its boundary, dc, is homeomorphic to a 2-sphere. 
The restriction of an SU(2) connection A to dc can be viewed as a connection 
A\g c on an SO (3) bundle over dc. There is then an SO (3) gauge transformation 
mapping (A, <&)|a c to (A, \^\as)\dc- Projecting A to a Cartan subalgebra of su(2) 
one obtaines a U(l) connection, a. This projection is called "abelian projection". 
It has been introduced by 't Hooft [36] to express Yang-Mills theories in terms 
of abelian gauge fields, charges and monopoles. We identify the magnetic charge 
contained inside c with the first Chern number of the curvature, F a , of a. 
The relation between $, A and a can be described in terms of the homotopy exact 
sequence 

~ n 2 (SO(3)) -> n 2 (S 2 ) -^7Ti(L/(l)) -^tti(50(3)) -> n^S 2 ) ~ 

where 

d : n G Z ~ 7r 2 (5 2 ) -> 2n G Z ~ tti(C/(1)) 

and 

i* : n E Z ~ 7Ti(t/(l)) -> nmod 2 G Z 2 ~ tti(50(3)) 

The group 7ri(C/(l)) classifies C/(l)-bundles over the (topological) 2-sphere S 2 , 
here identified with dc. The integer n G 7Ti(£/(l)) is given in terms of a connection 
a on a principal £7(l)-bundle over 9c by 




where p denotes faces of the cube c. 

The group ni(SO(3)) classifies S'0(3)-bundles over 5" 2 ~ dc. We have the fol- 
lowing relation between n mod 2 G ni(SO(3)) and a connection A on a principal 
SO (3) bundle over dc: 

e *™ = exp[z{arg TrP(e^„ A(P) )}], (5.5) 
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where is the connection 1-form of an SU{2) bundle over p obtained by lifting 
the 5'0(3)-bundle over p. If e t7Tn ^ 1 one cannot extend the £0(3)-bundle to 
the interior of the cube c. This signals the presence of an odd number of singular 
Z 2 -monopoles of SO (3) inside c. 

If $ is a map from S\ to S 2 of Kronecker index N then the first Chern num- 
ber corresponding to <9[<&] is 2N. Hence 't Hooft-Polyakov monopoles have even 
magnetic charge and, with an abuse of language, following [37], the term 't Hooft- 
Polyakov monopole stands for an arbitrary monopole configuration with even mag- 
netic charge. The most general configuration inside a cube c in R 3 is a combination 
of 't Hooft-Polyakov monopoles and Z 2 monopoles characterized by odd magnetic 
charge. 

Next we explain how to define the magnetic charge on the lattice by means of the 
abelian projection, following [38]. 

In model C, we start by choosing a unitary gauge by means of an SO(3) gauge 
transformation W. In this gauge the fields are given by 

= W^Wf 1 = (7 3 

9<ij> = Wig^yW' 1 (5.6) 
We write the S'L r (2)-valued gauge field g<ij> as a product 



9<ij> — C < ij > u < ij > (9) 



with 



C 



<ij> 



(1 - |c<ij>| 2 ) 2 ~ C <ij> 

C <ij> (1 ~~ \ C <ij> | 2 ) 2 

tl<y>(0) = , (5.7) 

where c < ij > G C, c< i:?> denotes its complex conjugate and e< ^ > = arg(^ <ij> )n. 
Under a t/(l)-gauge transformation {e~ lAj ^} of the original fields, 9 and c trans- 
form as 



9^9 + dA , c<ij> -> c <ij> e2 Ai e2 Aj 

Hence 9 is a £7(l)-gauge field, and c is a charged field of charge 1. 
We define the magnetic charge in a lattice cube c by 
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pEdc 

where dO is restricted to the range [— 2n, 27r]. Equation (5.8) is the lattice analogue 
of eq. (5.4). A Z 2 - charge in a cube c is defined by 

e i7rz c (g) = e * E peSc ^sx(ga P ) _ 

A plaquette p where 

e *argx(ffa p ) _ _2 

can be identified as the location of a Dirac string of a Z 2 -monopole intersecting 
the plane containing p. 

The relation established between (5.4) and (5.5) can be translated to the lattice 

as 

e inz c (g) _ e i-Km c (6) _ (5.10) 

With the help of the abelian projection, the phase transitions in model C have 
been analyzed numerically [39] in terms of condensations of magnetic currents. 

5.2 Monopole Green functions 

In model C, the magnetic charge of 't Hooft-Polyakov monopoles is even and 
the electric charge of the matter field is integer. Hence we are facing a situation 
analogous to the one discussed in sect 4.3, for q = 2, in models of class B. Following 
the ideas in that section, we propose to construct the two-point function of a 
monopole of magnetic charge 2 by summing over pairs of fluctuating strings of 
magnetic charge 1 with end points at the location of the monopole. 
As discussed in the previous section, these strings can be identified as Dirac strings 
of Z2-monopoles. They can be introduced by means of 't Hooft disorder fields. 
For a 2-surface E bounded by a lattice loop L, the disorder field is defined by 

D(E) = e -[So( 9 e^)-s ( 9 )] (5 _ n) 

where Sq is given in (2.8), E is a 2-form with values in {0, 7r}, whose support is 
dual to the surface E. 

The expectation value of -D(E) depends only on L. Clearly e lS<T3 takes values in 
{0, 1} ~ Z 2 , which is the center of the gauge group SU(2) . 
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In order to construct the monopole two-point function we start, by introducing a 
U(l) scalar field x on the sublattice A(i?, t). The field x is minimally coupled to a 
real gauge field A with coupling constant it. The action for x is given by S\(x, r), 
a functional defined in (4.25), with q = 2. 

Since «Sa(x> r ) is periodic in A with period 2, (e lXa; e _JXH +)A( fl)X o)(A) is the Fourier 
transform of a complex measure dfiR + (j x ) on the space of currents j x with values in 
7rZ which are constrained by 8j x = 2ti(5 x — 5 R+ ); see formula (4.26). We propose 
to define the two-point function for a 't Hooft-Polyakov monopole of magnetic 
charge 2 by 

G*{x,y) = lim c^c^i / d^ R+ (j x )dfi R _(j y )D LJ (j x , j y )) (5.12) 

R— >oo J 

with E constrained by 



*dT,(j x ,j y ;uj) = j x + j y + u + 2tt6 <r+r _ > , (5.13) 

where is a 3-form taking values in {0, 2tt} whose support is dual to a path 
connecting x to y. 

As the definition (5.12) involves a 't Hooft disorder field, the expectation value on 
the r.h.s. is independent of the choice of E, given u, and, since u takes value in 
{0, 27r} , also of u. 

In Appendix B it is shown, using the abelian projection, that the disorder field 
Du}(jx,j y ) introduces a source for 't Hooft-Polyakov monopole currents. For this 
purpose, a representation of < D u! {j x ,jy) > as a sum over configurations of mag- 
netic currents of even charge, with boundary given by {x, y}, is derived. This 
representation makes precise the idea that, in the abelian projection, the effect of 
the disorder field is to introduce a shift of d9 by 2E(j x , j y ; uj). Hence, according 
to definition (5.8), the magnetic charge in a cube c is changed by 2(<iE) c . 

Remark 5.1 In (5.10) one may replace the matrix 03 by $ x at the plaquette 
p = (x; hi, 112)1 since, in the abelian projection, $ x is replaced by 03, and the 
computation of the magnetic charge is then unchanged. 

Remark 5.2 The above construction cannot be adapted to obtain Green functions 
for the Z2-monopoles. For such Green functions, the current uj would take values 
in {0, 7r}, and the expectation value of the disorder field D^jxtjy) would depend 
on uj, hence on the Dirac strings of the Z2-monopoles, which is unphysical (see 
Appendix B). 



40 



One can argue that, in the Coulomb phase of model C correlation functions of non- 
zero total magnetic charge vanish (as proved at (3h = oo, where the model reduces 
to a U(l) theory), so that we expect the appearence of superselection sectors 
labelled by an even magnetic charge . According to the previous section, one may 
identify these sectors as 't Hooft-Polyakov monopole sectors. In the confinement 
phase of the model, Green functions of non-zero total charge are expected not to 
vanish due to "monopole condensation" (this is proved at (3h = oo). The large 
distance behaviour of the two-point monopole Green function could then be used 
to detect phase transitions in the Georgi-Glashow model, as in the U(l) gauge 
theory. 

5.3 Continuum 

In the final section we sketch how to define Green functions for quantum 't Hooft- 
Polyakov monopole fields in the formal continuum limit of the Georgi-Glashow 
model. 

It turns out that it is more convenient to work in a first-order formalism for the 
Yang-Mills field. (A formal relation of this formalism to ordinary Yang Mills and 
BF theories is discussed in [40]). 

We introduce a 2-form B with values in stt(2), which is a section of the vector 
bundle of 2-forms associated to Psu(2)- We replace the Yang-Mills euclidean 
action 

So(A) = ±J \F A \ 2 d 4 x 

by 

S (A, B) = i J Tr (Fa AB) + ^J \B\ 2 d 4 x (5.14) 

In the functional integral approach, integrating over B with the (white-noise) 
gaussian measure corresponding to the second term in (5.14), formally yields the 
standard weighting factor, e~ s °( A \ for Yang-Mills fields. 
In the first order formalism, we use 

DuUxJy) = [e 4 / MA * EW "* :w) U (5-15) 

as our 't Hooft disorder field, with the notations of the previous section adapted 
to the continuum. The notation [-] ren indicates that an ultraviolet multiplicative 
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renormalization is necessary. For the lattice theory, (5.15) corresponds to the 
choice mentioned in Remark 5.1. 

We wish to comment on the relation between our construction of monopole Green 
functions, based on expression (5.15) for the disorder operator, and the conven- 
tional semi-classical analysis of 't Hooft-Polyakov monopoles. For this purpose, 
we consider a monopole two-point function with y = rx (where, as usual, r denotes 
time reflection), and x = (x°, 0). We are interested in the asymptotic behaviour of 
this Green function, as x° becomes large. One strategy to analyze this behaviour is 
to attempt to evaluate the functional integral defining the Green function with the 
help of semi-classical techniques. In a semi-classical approximation, an evaluation 
of the expectation value 

lim ( / dfi R + (j x )dfj, R _(jy)D LJ (j x ,j y )) 

R— >oo J 

of the disorder operator in the formal functional measure of the Georgi-Glashow 
model is accomplished through an expansion of the functional measure around a 
back-ground field configuration. In the unitary gauge of the abelian projection, 
the disorder operator creates a "mean background gauge field" with curvature 
Fa(z) approximately given by 

a 3 B{z) 

for z = (z°,z), with \z°\ « x° : where B{z) is the rotation-covariant, static 
magnetic field generated by a magnetic monopole of magnetic charge 2, located at 
the origin. 

If, instead of the unitary gauge, we choose a gauge with the property that 

$( z ) = (5.16) 
\z\ 

the mean background gauge field has a field strength F^(z) approximately given 
by 

B H ~ P (z) = (5.17) 
for z = (z ,z), with \z \ << x , \z\ large. 

The field in (5.17) describes the large-distance behaviour of the field strength of 
the classical 't Hooft-Polyakov monopole solution [5]. 
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We should comment on the notion of "mean background (gauge) field" used in the 
arguments above: Every field configuration contributing to the expectation value 
of a disorder field Duijxijy) actually satisfies the non-abelian version of Dirac's 
quantization condition [28,37]. Such a field configuration is therefore singular near 
the support of the magnetic flux currents, j x and j y . However, after integrating 
over j x and j y with the complex measure dnn + (j x )dfiR_(j y ) and taking the limit 
R — > oo, and after integrating out the high-frequency (short-distance) modes of 
the fields, i.e., after ultraviolet renormalization ("coarse graining"), the resulting 
background field configuration approaches the one described above. This is be- 
cause the constraint of flux quantization of the Dirac strings is softened under 
renormalization and is actually expected to scale to zero in the limit of very large 
distances scales. This phenomenon is analogous to the one discussed in the abelian 
models of class B. 

Taking it for granted that, on large distance scales, the background field config- 
uration described in (5.16), (5.17) dominates the functional integral appearing in 
the (numerator of the) monopole Green function G(x,rx), as x° becomes large, 
the relation of our approach to the conventional semi-classical analysis of (e.g., 
the mass of) the 't Hooft-Polyakov monopole becomes clear. 
We conclude this section with a comment on dyons. If a topological term 



is added to the action of the Georgi- Glashow model then the excitation spectrum 
of the model contains dyons . 

Dyon Green functions can be constructed with the help of disorder fields. The con- 
struction is analogous to the one discussed in the context of the Cardy-Rabinovici 
model, in sect. 4.5. If is not an integer multiple of 2tt dyons carry a fractional 
charge. As approaches an integer multiple of 2n, the dyons correspond to those 
first described by Julia and Zee [41]. The spectrum of the theory is periodic in 0, 
with period 2tv. 

Appendix A (Gaussian evaluation of F(E\dA)) 

[In this appendix all symbols referring explicitly to the sublattice A(R, t) are omit- 
ted, e.g., we write d instead of d\ etc..]. 

To analyse F(E\dA) it is convenient to work on the dual lattice and introduce an 
auxiliary gauge field, C: 
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/ rw dc <ij>e - wii^n Ep;fp=0 e -i(c,-^) e i (d c-- g , g (p)) 



where, comparing with (4.18), p =* v and u(p) is a fixed, integer-valued solution 
of 

Sa(p) = p. 

Following the techniques developed in [12] (see also [2,3]), we divide a configuration 
of currents, p, into (not necessarily connected) networks, {p a }, such that the 
distance d(p a , pp) > 2, a ^ (3. 

We now renormalize the activities of the current networks, {p Q }, using the method 
of complex translations. 

Let B(p a ) be a set of links in the support of p a such that two links in B{p a ) do 
not belong to a common plaquette and such that 

\Pce<ij>\ 2 > c {Pai Pa) 5 

<ij>eB( Pa ) 

where 

c" 1 = card {£' : £' ^ £, {£, £'} G dp for some p} = 18. 
We set 



P' a = Pa\t3( Pa ) p' = ^2p'c 

and perform the complex translation 



a 



C + i^p' + i(3 x A~ l *dA 
n x 



where 



n = card{p : dp G< ij >} = 6 
As a result every current network p a obtains a complex activity 

Z A {p a , C, E) = e -^l\K\\ 2 +2(P^*dA) ]el ( d C-*E,a( P) ) 

in the numerator of (A.l), and an activity ZA(pa] C) = ZA(p a , C, 0) in the denom- 
inator of (A.l), where C = C- ^{5dC)\ p > a . 
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If dA is sufficiently small so that 



\z A ( Pa ,C)\,\z A (Pa,C,E)\ < e-^U?^ 2 (A.2) 

the current networks form a dilute gas, and we can exponentiate their contributions 
in the form of a Mayer series: 

F(E\dA) = J n<lJ> — = ^ (A.3) 

/n< <i >dc< y >e-^ l|dcr||a+ ^ or( ^ (£ » c) 

where £ is a collection of current networks p in which a single p can occur an 
arbitrary number of times , aj , (£) is a standard combinatorial factor which enforces 
all C to have connected supports and 

z A (C,C,E)= l[z A (p,C,E). 

pec 

We give an approximate evaluation of (A.3) by keeping only the leading contribu- 
tion in C, which correspond to plaquette terms, and expanding the activities z A 
to second order in C, A, E. 
F(E\dA) is then approximated by 



rrr rir ^-wll^H^E, e-°^Hi + ^{*A T )i)(dc-*E)i +l p x (*A T ) p {dc) p }} 
J ll<ij> a<^<ij>e x * 

rrr Ar S-^W dC \V-Y. v e- ^x)(l+^(* A T)2) {d c)2 + if3 x ^ A T )p{d c) p ]} 

J ll<ij> a^<ij>e Px p 

(A.4) 

where A T = 5A _1 cL4, and we have used that 

(dC*E) = 0. 

Explicit evaluation of (A.4) to second order in A gives 

F G (E\dA) = exp[-e-°^)(-||£;|| 2 + ^\\A T ■ E\\ 2 )] 

Appendix B 

In this appendix we exhibit a representation of (D u (j x , j y )) in model C in terms 
of magnetic currents. We start by replacing the SU (2) gauge field g with a couple 
of new variables {U, a}. U is the gauge coset variable given on a link < ij > by 
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(7<ij>r, where V ~ Z 2 is the center of the gauge group SU(2); a is a 2-form with 
values in {0, n} ~ Z 2 . 

The variables t/ and cr are not completely independent: it can be proved [42] that 
e i-Kz c { g ) d e f] nec i j n (5_g) j s a function of the coset field U, which we denote by 
e iirz c (u)^ anc j following constraint holds: 

e iirz c (U) _ e in(da) c (B.l) 

for every cube c. 

The partition function of model C can be rewritten as 

z =nj n ^<«>n d *< 

o <ij> i 



( J2 p M(u 9p )e^-s 1 (u,'P) Y[s(e i[7VZc{u) - {da)c] ) (B.2) 

c 

where 



|X|(^ P ) = |XW|. (B.3) 

The introduction of the disorder field -D(S) induces a shift of a v by E p . We perform 
a duality transformation in the Z 2 variable. Let r denote the 1-form with values 
in {0, 1} dual to cr. Defining 

S (U;dr) = -J]lncosh[| X |(t/a p )] +* (dr) p \ntanh[\ X \(U dp )}, 
p 

the expectation value of 't Hooft disorder field can be written 



T <ij> i c 

(B.4) 

To exhibit the wordlines of magnetic currents we perform the abelian projection. 
We can decompose U <i j > as in the l.h.s. of (5.7), but with 6 < ij> restricted to the 
range (— tt, tt). 

As a consequence of £7(l)-gauge invariance, we can expand 

/ II rfC <lJ>e -[ Sl ( c ^) +So ( CeICT3f; ^ (B.5) 

J <ij> 
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as a Fourier series in dd. The Fourier coefficients are denoted by F(n; dr), where 
n is an integer valued 2-form. 
Defined the 1-form £ by 



5n = £ (B.6) 

we decompose the 2-form n as 

n = n[£] +* d$, 

where n[£] is an integer-valued solution of (B.6) and £ a Z-valued 1-form in the 
dual lattice. Furthermore, we define a Z/2-valued 1-form a in the dual lattice by 

a = i+\r (B.7) 

and we adopt the notation: 

F(n[£] + X; dr) = F(n[£]\*da) (B.8) 
Making use of eq.(5.10), we obtain 



= E F(n[£]\*da) 



I 



Z 

[a] e-.se=o 

H de^ye^^+^^e^ (B.9) 

<ij> 



where [a] denotes a gauge equivalence class of a. 

In eq.(B.9) we can replace a by a real-valued 1-form A by inserting the term 
Y^p-8 p =o e l ( A,p \ where p is a 4n Z-valued 1- form. 

We split m(6>) into a component of even magnetic charge m e {9) and a component 
of odd magnetic charge m (9). Shifting p by 2itm e {Q), we obtain the following 
identity: 

(D»U*Jv)) = ^ E / d M E F(n[t]\*dA) I Hd9 <lJ>e ^ 

£:5e=0^ p:8p=0 J <ij> 

e i(A,2j x +2j y +2aj+p+2nm (9)+4:TTS <R+R _ > ) ^Q^) 

where d[A] denotes formal integration over gauge equivalence classes of A . In 
eq.(B.lO) worldlines of 't Hooft-Polyakov monopoles are described by 2uj + p and 
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they exhibit sources at {x} and {y}. Analogously 2nm describe the virtual tra- 
jectories of Z 2 monopoles. The representation (B.10) shows also explicitely the 
independence of the choice of to in the construction of Green functions of 't Hooft- 
Polyakov monopoles. 

Taking u> with values in {0, 7r}, a similar representation shows that our construction 
of Green functions cannot be adapted to Z2-monopoles, because a change in u> 
cannot be reabsorbed by field redefinition. 

Acknowledgements 

One of us (P.M.) gratefully acknowledges A. Di Giacomo and D. Marenduzzo for 
useful discussions and E.T.H. in Zurich for kind hospitality. 

This work was partially supported by the European Commission TMR programme 
ERBFMPX-CT96-0045 to which P.M. is associated. 

References 

[1] P.A.M. Dirac, Can. J. Phys. 33, 650 (1955). 

[2] J. Frohlich, P.A. Marchetti in "The Algebraic Theory of Superselection Sectors. 

Introduction and Recent Results". D. Kastler ed., World Scientific 1990. 
[3] J. Frohlich, P.A. Marchetti, Europhys. Lett. 2, 933 (1986) and unpublished notes. 
[4] N. Seiberg, E. Witten , Nucl. Phys. B426, 19; B431, 484 (1994). 
[5] G. 't Hooft, Nucl. Phys. B79, 276 (1984); A.M. Polyakov JETP Lett. 20, 194 

(1974). 

[6] B. Julia, A. Zee, Phys. Rev. Dll, 2227 (1975). 
[7] A.M. Polyakov,Nucl. Phys.B120; 429 (1977). 

[8] H.A. Kramers, G.H. Wannier, Phys. Rev. 60, 252 (1941); F. Wegner, J. Math. 

Phys. 12, 2259 (1971); G. 't Hooft, Nucl. Phys. B138, 1 (1978). 
[9] S. Mandelstam, Ann. Phys. 19, 1 (1962). 
[10] G. 't Hooft, Nucl. Phys. B138, 1 (1978); Nucl. Phys. B153, 141 (1979). 
[11] J. L. Cardy, E. Rabinovici, Nucl. Phys. B205, 1 (1982); J.L. Cardy, Nucl. Phys. 

B205, 17 (1982). 
[12] J. Frohlich, T. Spencer, Commun Math. Phys. 83, 411 (1982). 
[13] A. Guth, Phys. Rev. D21, 2291 (1980). 
[14] D. Brydges, E. Seiler, Journal Stat. Phys. 42, 405 (1986). 

[15] E. Seiler, "Gauge Theories as a Problem in Constructive Quantum Field Theory 
and Statistical Mechanics", Lecture Notes in Physics 159, (1982). 



48 



[16] see e.g. C. Borgs, F. Nill, "The Phase Diagram of Abelian Lattice Higgs Model; 

A Review of Rigorous Results". Preprint ETH 87-0069; J. Bricmont, J. Frohlich, 

Nucl. Phys. B230, 407 (1984). 
[17] J. Frohlich, P.A. Marchetti, Commun. Math. Phys. 121, 177 (1989); V. Miiller 

Z. Phys. C51, 665 (1991). 
[18] R. Brower et al., Phys. Rev. D25, 3319 (1982). 

[19] J. Frohlich, P.A. Marchetti, Commun. Math. Phys. 112, 343 (1987). 
[20] J. Frohlich, P.A. Marchetti, Commun. Math. Phys. 116, 127 (1988). 
[21] J. Frohlich, P.A. Marchetti, Commun. Math. Phys. 121, 177 (1989). 
[22] E.C. Marino, B. Schroer, J.A. Swieca, Nucl. Phys. B200, [FS 14], 472, (1982); K. 

Fredenhaegen, M. Marcu, Commun. Math. Phys. 92, 81 (1983); K. Szlachanyi, 

Commun. Math. Phys. 108, 319 (1987); J. Frohlich, P.A. Marchetti, Lett. Math. 

Phys. 16, 347 (1988); Nucl. Phys. B355, 1 (1990); P.A. Marchetti, Europhys. 

Lett. 4, 633 (1987). 

[23] J. Bricmont, J. Frohlich, Nucl. Phys. B251, [FS 13], 517 (1985); Commun. Math. 

Phys. 98, 553 (1985); P.A. Marchetti, Commun. Math. Phys. 117, 501 (1988). 
[24] J. Frohlich in "Scaling and Self-similarity in Physics", J. Frohlich ed., "Progress 

in Physics" , Birkhauser, 1983. 
[25] T. Kennedy, C.King, Phys. Rev. Lett. 55, 776 (1985); Commun. Math. Phys. 

104, 345 (1986). 
[26] D. Buchholz, Commun. Math. Phys. 85, 49 (1982). 

[27] L. Polley, U.J. Wiese, Nucl. Phys.B356, 629 (1991); M.I. Polikarpov, L. Polley, 
U.J. Wiese, Phys. Lett. B253, 212 (1991); L. Del Debbio, A. di Giacomo, G. 
Paffuti, Phys. Lett. B349, 513 (1995). 

[28] S. Coleman in "The Unity of the Fundamental Interactions" (Erice 1981), A. 
Zichichi ed., Plenum Press, 1983. 

[29] D. Durhuus, J. Frohlich, Commun. Math. Phys. 75, 103 (1988). 

[30] D. Marenduzzo, Tesi di laurea, University of Padova. 

[31] J. Frohlich in "Unified Theories of Elementary Particles. Critical Assessments and 
Prospects". P. Breitenhohner, H.P. Durr eds., Springer Lecture Notes in Physics 
vol. 160 (1982). 

[32] D. Brydges, J. Frohlich, T. Spencer, Commun. Math. Phys. 83, 123 (1982). 

[33] L. Del Debbio, A. Di Giacomo, G. Paffutti, L. Pieri, Phys. Lett. B355, 255 
(1995); A. Di Giacomo, G. Paffuti, Phys. Rev. D 56, 6816 (1997); A. I. Veselov, 
M.I. Polikarpov, M.N. Chernodub, JETP Letters 63, 411 (1996). 

[34] E. Witten, Phys. Lett. B86, 283 (1979). 



49 



[35] A. Jaffe, C. Taubes "Vortices and Monopoles. Structure of a static gauge theories" . 

Progress in Phys. 2; A. Jaffe, D. Ruelle eds., Birkhauser 1980. 
[36] G. 't Hooft, Nucl. Phys. B120, 455 (1981). 

[37] F. Englert in "Hadron Structure and Lepton-Hadron Interactions" (Cargese 1977) 

M. Levy et al. eds., Plenum Press, 1979. 
[38] A.S. Kronfeld, G. Schierholz, U.J. Wiese, Nucl. Phys. B293, 461 (1987). 
[39] A.S. Kronfeld, M.L. Laursen, G. Schierholz, U.J. Wiese, Phys. Lett. B198, 516 

(1987). 

[40] A.S. Cattaneo et al., "Four-Dimensional Yang-Mills theory as a Deformation of 

Topological BF theory" ; [hep-th/9705123| , to appear in Commun. Math. Phys. 
[41] B. Julia, A. Zee, Phys. Rev. Dll, 2227 (1975). 
[42] G. Mack, V.B. Petkova, Z. Phys. C12, 177 (1982). 



50 



Figures 

(One space dimension is suppressed) 




Fig. 1. 
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Fig. 3. 
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